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The possibility that the observed cosmological baryon density might 
be a consequence of non-perturbative dynamics at the electroweak phase 
transition is one of the most exciting topics at the interface of particle 
physics and cosmology. However, determining the viability of electroweak 
baryogenesis in (minimal extensions of) the standard model requires knowl- 
edge of the equilibrium behavior of the electroweak phase transition, non- 
equilibrium dynamics around expanding bubble walls, and non-perturbative 
baryon violating processes in both the high and low temperature phases [1]. 

This talk will examine the current understanding of the electroweak 
phase transition. Determining the order of the phase transition, and the 
magnitudes of the latent heat, correlation lengths, and the baryon violation 
rate at the transition are some of the key questions. Viable baryogenesis 
scenarios require a first order transition with rapid suppression of baryon 
violating transitions in the low temperature phase [1, 2, 3]. 

Performing reliable quantitative calculations of electroweak transition 
properties is challenging, in part because there is important dynamics, both 
perturbative and non-perturbative, on many diff^ering length scales. Ap- 
proximation methods which have been applied to this problem include weak 
coupling perturbation theory (or mean field theory) [2, 4, 5], e-expansions 
[6, 7, 8], and numerical simulations [9, 10, 11, 12]. Each of these approaches 
have significant (but differing) limitations; for example, ordinary perturba- 
tion theory is valid if the Higgs mass is much lighter than the VF-boson 
mass, but does not appear to be trustworthy for the experimentally al- 
lowed range of possible Higgs masses. Nevertheless, all of these approaches 
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have helped to elucidate the rich interplay of the physics at the electroweak 
phase transition. 

I will focus on the behavior of the minimal standard model, despite the 
fact that electroweak baryogenesis in the strictly minimal model cannot ex- 
plain the observed cosmological baryon density. There are two basic prob- 
lems with baryogenesis in the minimal model. For experimentally allowed 
values of the Higgs mass, most of the baryon excess which is produced at the 
phase transition is destroyed shortly thereafter [2]. And even if that weren't 
the case, the baryon excess produced in the minimal standard model would 
still be far too small to agree with present day observations [1]. The first 
problem will be discussed further below. The second problem is a conse- 
quence of the nature of CP violation in the minimal model. Without funda- 
mental CP violation, no baryon excess can be dynamically produced. But 
in the minimal standard model, CP violation arises only through a phase 
in the CKM matrix elements, and can only appear in processes involving 
all three generations of fermions. Because of the freedom to redefine the 
relative phases of fermion fields, one may show that any CP violating tran- 
sition must involve a reparamctcrization invariant combination of Yukawa 
couplings and mixing angles which is tiny [3], 



Given this suppression, simple estimates show that it is essentially impossi- 
ble to produce a baryon-to-photon ratio of the required 10"^'^ magnitude [1]. 

Both of these difficulties can be avoided in extensions of the minimal 
standard model which incorporate, for example, a single additional scalar 
field. This allows explicit CP violation in the scalar sector, and can increase 
the sphaleron mass and thereby slow down the rate of baryon- violating re- 
actions just after the transition. However, the challenges involved in analyz- 
ing the electroweak phase transition are, for the most part, generic to any 
electroweak theory. Consequently, the minimal standard model provides 
a useful toy model (with the fewest adjustable parameters) which may be 
used as a testing ground for quantitative calculations of electroweak phase 
transition properties. 

1. Perturbation Theory 

Since the standard model is weakly coupled (at scales of several hundred 
GeV), perturbation theory is the most obvious approach for studying the 
electroweak phase transition. A one-loop calculation of the effective poten- 
tial for the Higgs field is elementary and gives (schematically) 



^CP < 10 
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= (-//^ + a g^T^) (f^ -h g^T 4>^ + \4>* + (^-independent) (3) 

For simplicity, only the thermal VF-boson contribution is indicated in (2). 
The VF-boson energy E{k) = y/F^TMw(W with Mw(0) = \g4>, and a 
and h are dimensionless constants which I will henceforth suppress. 




Figure 1. The form of the free energy, as a function of ^, for different temperatures. 

The 0{T'^) thermal correction to the quadratic term in the potential is 
responsible for driving the transition from the Higgs phase, with {(f>) ^ 0, to 
an unbroken symmetry phase with {(f)) = 0, as the temperature is increased. 
If corrections to this one-loop result are negligible, then the presence of 
the 0{g^T) cubic term will cause the transition to be first-order, with a 
barrier separating co-existing broken and unbroken symmetry phases at 
the transition temperature T^. (Sec Fig. 1.) In the Higgs phase near the 
transition, the quadratic, cubic, and quartic terms in the potential (3) are 
all comparable in size. Consequently, (f)c ~ g^T/X, or 

(4) 

T=0 

These ratios are important for several reasons. First, the rate of non- 
perturbative baryon violating reactions in the low temperature Higgs phase 
is controlled by a thermal activation energy given by the (temperature 
dependent) sphaleron mass, 

rAB~T4e-^=p>>(^V^, (5) 

and the sphaleron mass Msph(r) equals IGir M^^(T)/g^ (to within a factor 
of 2). Hence, the relation (4) implies that the baryon violation rate, just 
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after the transition, is exponentially sensitive to the zero temperature ratio 
of Higgs and W^-boson masses, 



InFAB--^^ 



Ml 



(6) 

r=o 



In order for electroweak baryogenesis to be viable, baryon violating reac- 
tions after the transition must turn off sufficiently rapidly so that the baryon 
asymmetry produced during the transition can survive to the present day. 
The simple result (6) implies that this will only be possible if the Higgs 
to W mass ratio is sufficiently small. More detailed analysis, based on the 
one-loop effective potential in the minimal standard model, produces an es- 
timate of 35 Gev for the upper limit on the Higgs mass if baryon violating 
rates after the transition arc to be acceptably small [2]. This, of course, is 
inconsistent with the current experimental lower bound on the Higgs mass 
of about 65 GeV. 

As noted above, one can circumvent this "no-go" result if minor ex- 
tensions (such as adding a second scalar field) are made to the minimal 
standard model. However, one should first ask whether the 35 GeV bound 
in the minimal model, or any other conclusion based on the one-loop analy- 
sis (including the predicted first order nature of the transition!) is reliable. 
To answer this, one must understand whether higher order corrections to 
the one loop results are significant. 

Life would be simple if the temperature T were the only relevant scale 
for physics at the phase transition. If this were the case, then higher order 
corrections would automatically be suppressed by powers of (or aw) at 
the scale T, which is, in fact, small.^ However, life is not so simple. Con- 
sider, for example, any effective potential diagram containing only gauge 
field lines. The contribution of each loop may be estimated as^ 

g'^T I d\ (q2 + Mw(T)2) ~ g'^T/M^iT) . (7) 

Therefore the real loop expansion parameter is not g"^ , but rather equals 
g'^T /M^^{T). Consequently, the reliability of perturbation theory, in the 
Higgs phase near the transition, is controlled by the the same ratio of 



9"% Ml 



w 



T=0 



^This assumes that the Higgs is Ught enough to be weakly coupled. 

^Recall that Euchdean space frequencies are discrete at finite temperature, g° = Wn = 
(27rnT). Convergent diagrams are dominated by the static n=0 frequency component. 
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which governs the baryon violation rate. Perturbation theory is not rehable 
unless the physical Higgs is sufficiently light. 

Of course, the above estimate docs not determine whether perturba- 
tion theory really breaks down at Mh = Mw/2, or at Mh = 2Afw. For 
a quantitative answer, one must actually compute higher order corrections 
to various quantities of interest. Fortunately, a number of two-loop cal- 
culations of phase transition properties have now been performed. Fig. 2 
illustrates the result for the effective potential at a Higgs mass of 35 GeV 
[4, 5]. The magnitude of {4>) in the Higgs phase at the transition shifts by 
only about 20%, but the height of the barrier changes by almost a factor 
of three! 




20 40 60 80 100 120 140 
(t)[GeV] 



Figure 2. The effective potential at the critical temperature for Afn(O) = 35 GeV 
and mt(0) = 110 GeV. The dashed and solid lines are the one-loop and two- loop results 
respectively. [Why a 110 GeV top mass? Because this is an old graph. But the results 
aren't particularly sensitive to mt.] 

Since both the barrier height, and the Higgs expectation value, are gauge 
dependent, one might legitimately wonder whether truly physical quantities 
will be better behaved. However, a two-loop evaluation of the dimensionless 
latent heat, AQ/T^, and surface tension, cr/T^, at a Higgs mass of 50 GeV 
finds 50-100% changes over the one-loop results [12]. Consequently, there 
appears to be good reason to expect a Higgs mass somewhere between 35 
and 50 GeV to be the upper limit for the reliability of perturbation theory 
for most quantities.^ 

^Note, however, that the authors of ref. [11] claim that using the 3d renormalization 
group to improve some logarithmic corrections can delay the breakdown of perturbation 
theory. 
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Finally, even when the Higgs to W mass ratio is small, so that pertur- 
bation theory in the Higgs phase is well behaved, the reliability of pertur- 
bativc calculations of phase transition properties is limited by the presence 
of non-perturbative physics in the unbroken symmetry phase. The long 
distance physics of the unbroken phase is described by a three dimensional 
non-Abclian theory with a confinement scale of g^T, and has a free energy 
density of order {g^T)"' which is incalculable in perturbation theory. Con- 
sequently, for small values of the Higgs fields, where Mw(^) ^ 9^T, any 
perturbative approximation to the effective potential will have a system- 
atic uncertainty of order g^T'^. This will generate an uncertainty in the 
value of the transition temperature, or any other phase transition property. 
However, when X/g"^ is small, this uncertainty is of the same order as the 
unknown four- loop contributions to the free energy in the Higgs phase [4]. 

2. Numerical Simulations 

The electroweak phase transition can be also be studied numerically by 
performing stochastic simulations in a lattice version of the theory. Sim- 
ulating the fuh theory (with chiral fermions and SU(3) x SU(2) x U(l) 
gauge fields) is not practical. However, since one is interested in high tem- 
perature physics, all non-zero frequency Fourier components of fields have 
very short correlation lengths (of order T~^) and may be integrated out 
perturbatively. Fermion fields, being antiperiodic, have no static compo- 
nents and may be completely eliminated. Once the fermions are gone, the 
SU(3) gauge field is decoupled and may be dropped, as may the U(l) gauge 
field if the small weak mixing angle is treated as perturbation. This leaves 
an SU(2)-Higgs theory, which is straightforward to define on a lattice. 

During the past two years, several groups have performed large scale 
simulations of the finite temperature phase transition in SU(2)-Higgs the- 
ory. The authors of refs. [9, 12] have chosen to simulate a four-dimensional 
SU(2)-Higgs theory with a periodic "time" dimension consisting of Nt = 
2, 3, . . . lattice spacings. In contrast, the authors in refs. [10, 11] begin with 
the effective three dimensional theory which results from integrating out all 
non-static Fourier components. Spatial lattice sizes in these simulations are 
substantial (typically 16'^ to 32'''). 

I will not attempt to summarize all the results of these lattice simula- 
tions (which include data on the transition location and character, latent 
heat, surface tension, and correlation lengths). Interested readers should re- 
fer to the latest papers. However, if one asks what lessons can be extracted 
from the results of these simulations, the the following points appear to be 
fairly well established: 

a. Reasonably good calculations can be performed, at least for light Higgs 
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masses. Finite lattice size and non-zero lattice spacing errors are (at 
least beginning to be) under reasonable control. 

b. For Higgs masses below Mw, a first order phase transition is seen. The 
discontinuity at the transition becomes steadily weaker with increasing 
Higgs mass. 

c. When Mh < 20 GeV, the lattice results agree quite well with pertur- 
bative predictions. When Mh > 50 GeV, the strength of the transition 
appears to be larger than two-loop perturbative predictions. 

d. Simulations (of sufficient accuracy to study a weakly discontinuous 
transition) become progressively more difficult as the Higgs mass grows. 
So far, no convincing conclusions can be drawn about the nature of 
the transition when Mh > Mw- 



3. e-expansions 

The e-expansion provides an alternative systematic approach for comput- 
ing the effects of (near)-critical fluctuations. It is based on the idea that 
instead of trying to solve a theory directly in three spatial dimensions, it 
can be useful to generalize the theory from 3 to 4— e spatial dimensions, 
solve the theory near four dimensions (when e ^ 1), and then extrapolate 
to the physical case of 3 spatial dimensions. Specifically, one expands phys- 
ical quantities in powers of e and then evaluates the resulting (truncated) 
series at e = 1 [13]. This can provide a useful approximation when the 
relevant long distance fiuctuations are weakly coupled near 4 dimensions, 
but become sufficiently strongly coupled that the loop expansion parameter 
is no longer small in three dimensions. 

Scalar (p^ theory (or the Ising model) is a classic example. In four 
dimensions, the long distance structure of a quartic scalar field theory is 
trivial; this is reflected in the fact that the renormalization group equation 
li{d\/dii) = + O(A^) , has a single flxed point at A = 0. In 4— e dimen- 
sions, the canonical dimension of the field changes and the renormalization 
group equation acquires a linear term, 

^^ = _eA + cA2 + 0(A3). 
djJL 

This has a non-trivial fixed point (to which the theory flows as ^ decreases) 
at A* = e/c + 0(e^)- The fixed point coupling is 0(e) and thus small near 
four dimensions, but grows with decreasing dimension and becomes order 
one when e = 1. Near four dimensions, a perturbative calculation in powers 
of A is reliable and directly generates an expansion in powers of e. 

The existence of an infrared-stable fixed point indicates the presence 
of a continuous phase transition as the bare parameters of the theory are 
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varied. Performing conventional (dimensionally regularized) perturbative 
calculations and evaluating the resulting series at the fixed point, one finds, 
for example, that the susceptibility exponent (equivalent to the anomalous 
dimension of (jP') has the expansion [13, 14] 

7 = 1 + 0.167 e + 0.077 - 0.049 + O(e^) . (8) 

Adding the first three non-trivial terms in this series, and evaluating at 
e = 1, yields a prediction which agrees with the best available result to 
within a few percent [15, 16]. 

Inevitably, perturbative expansions in powers of A are only asymptotic; 
coefficients grow like n!, so that succeeding terms in the scries begin growing 
in magnitude when n > 0(1/A). Expansions in e are therefore also asymp- 
totic, with terms growing in magnitude beyond some order n > 0(l/e). If 
one is lucky, as is the case in the pure scalar theory, 0(l/e) really means 
something like three or four when e = 1 and the first few terms of the series 
will be useful. If one is unlucky, no terms in the expansion will be useful. 
Whether or not one will be lucky cannot be determined in advance of an 
actual calculation. 

To apply the e-expansion to electroweak theory, one begins with the 
full 3+1 dimensional finite temperature Euclidean quantum field theory 
(in which one dimension is periodic with period (5 = 1/T) and integrates 
out all non-static Fourier components of the fields. The integration over 
modes with momenta of order T or larger may be reliably performed using 
standard perturbation theory in the weakly-coupled electroweak theory. 
This reduces the theory to an effective 3-dimensional SU(2)-Higgs theory.^ 
The effective theory depends on three relevant renormalized parameters: 

51 (T)^ — the SU(2) gauge coupling, 
Ai(r) — the quartic Higgs coupling, 

mi(T)^ — the Higgs mass (squared). 

Next, one replaces the 3-dimensional theory by the corresponding 4— e 

dimensional theory (and scales the couplings so that gf/e and Ai/e are 
held fixed). This is the starting point for the e-expansion. When e is small, 
one may reliably compute the renormalization group flow of the effective 

^Fermions, having no static Fourier components, are completely eliminated in the 
effective theory. For simplicity, the effects of a non-zero weak mixing angle and the 
resulting perturbations due to the U(l) gauge field are ignored. Finally, one may also 
integrate out the static part of the time component of the gauge field, since this field 
acquires an 0{gT) Debye-screening mass. 



Figure 3. The renormalization group flow for an SU(2)-Higgs tlieory. Arrows indicate 
the direction of decreasing renormahzation point. The dashed hne is the trajectory which 
flows from an initial set of couplings {gf, Ai) into the region where X <^ g^. 



couplings. The renormalization group equations have the form 




(10) 



(9) 



The precise values of the coefficients (and the next order terms) may be 
found in reference [6] . These equations may be integrated analytically, and 
produce the flow illustrated in flgure 3. 

Note that a non-zero gauge coupling renders the Ising fixed point at A = 
e/c unstable, and that no other (weakly coupled) stable renormalization 
group fixed point exists. Trajectories with > eventually cross the 
A = axis and flow into the region where the theory (classically) would 
appear to be unstable. Such behavior is typically indicative of a first-order 
phase transition [17]. To determine whether this is really the case, one 
must be able to perform a reliable calculation of the effective potential 
(or other physical observables) . As discussed earlier, the loop expansion 
parameter for long distance physics is \{p) / g'^{ii). Consequently, the best 
strategy is to use the renormalization group to flow from the original theory 
at fj, = T, which may have \{T)/g'^{T) large, to an equivalent theory with 
jjL for which Xi^jj) / g^ {jj) is small. This is equivalent to the condition 
that one decrease the renormalization point until it is comparable to the 
relevant scale for long distance physics, specifically, the gauge boson mass, 
M. By doing so, one eliminates large factors of [{M/jiY — l]/e which 
would otherwise spoil the reliability of the loop expansion. (This, of course. 
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is nothing other than the transcription to 4— e dimensions of the usual 
story in 4 dimensions, where appropriate use of the renormahzation group 
allows one to sum up large logarithms which would otherwise spoil the 
perturbation expansion.) 

For small e, the change of scale required to flow from an initial the- 
ory where Xi/gl = 0{1) to an equivalent theory with X{iJ,)/g^{^) ^ 1 is 
exponentially large; the ratio of scales is 

Given the parameters g^{n), A(/x) and m?{fi) of the resulting effective 
theory, one may use the usual loop expansion to compute interesting physi- 
cal quantities. Because the change in scale is exponentially sensitive to 1/e, 
the result for a typical physical quantity will have the schematic form 

= f[9\,,),X{f^),m\f,)](^^Y (11) 



e* (1 + 0(e) + • • •) exp 



^ (1 -I- 0(e) -t- 



(12) 



In general, a calculation accurate to 0{e"') requires an n-loop calculation 
in the final effective theory, together with n+1 loop renormahzation group 
evolution. 

To obtain predictions for the original theory in three spatial dimensions, 
one finally truncates the expansions at a given order and then extrapolates 
from e <C 1 to e = 1. Just as for the simple (f)'^ theory, the reliability of the 
resulting predictions at e = 1 can only be tested a-posteriori. 

This procedure has been carried out for a variety of observables char- 
acterizing the electroweak phase transition at both leading and next-to- 
leading order in the e-expansion [6]. Leading order results are available for 
the scalar correlation length in both symmetric and asymmetric phases, 
the free energy difference AF{T) between the symmetric and asymmetric 
phases, the latent heat AQ = —T{dAF/dT)\T^, the surface tension a be- 
tween symmetric and asymmetric phases at Tc, the bubble nucleation rate 
rN(T) below Tc, and the baryon violation (or sphalcron) rate rAB(7c). 

The lowest order e-expansion predictions differ from the results of stan- 
dard one-loop perturbation theory (performed directly in three space di- 
mensions) in several interesting ways. First, the e-expansion predicts a 
stronger first order transition than does one-loop perturbation theory (as 
long as Mh < 130 GeV). The correlation length at the transition is smaller, 
and the latent heat larger, than the perturbation theory results. The size 
of the difference depends on the value of the Higgs mass (see reference [6] 
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for quantitative results). Naively, one would expect that a stronger first 
order transition would imply a smaller baryon violation rate, since a larger 
effective potential barrier between the co-existing phases should decrease 
the likelyhood of thermally-activated transitions across the barrier. This 
expectation is wrong (in essence, because it unjustifiably assumes that the 
shape of the barrier remains unchanged) . Along with predicting a strength- 
ing of the transition, the e-expansion predicts a larger baryon violation rate. 
This occurs because the baryon violation rate is exponentially sensitive to 
the sphaleron action (or mass). 

Tab OC exp — S'sphaleron , 

and the sphaleron action depends inversely on e, 

'S'sphaleron ~^JJx) ^^^1^^ ' 

Hence, unlike other observables, the exponential sensitivity to 1/e in the 
baryon violation rate does not arise solely from an overall power of the scale 
factor /x/T. 

Note that an increase in the baryon violation rate (compared to standard 
perturbation theory) makes the constraints for viable electroweak baryo- 
gencsis more stringent; specifically, the (lowest order) e-expansion suggests 
that the minimal standard model bound Mh < 35-40 GeV derived using 
one-loop perturbation theory in ref. [2] should be even lower, further ruling 
out electroweak baryogenesis in the minimal model. 

As emphasized earlier, in general there is no way to know, in advance 
of an actual calculation, how many terms (if any) in an e-expansion will 
be useful when results are extrapolated to e = 1. Therefore, to assess 
the reliability of an e-expansion one must be able to test predictions for 
actual physical quantities. For the electroweak theory, two types of tests 
are possible: 

A. X g^. In the limit of a light (zero temperature) Higgs mass, or 
equivalently small Xi/gf, the loop expansion in three dimensions is re- 
liable. Hence, although this is not a realistic domain, one may easily 
test the reliability of the e-expansion in this regime by comparing with 
direct three-dimensional perturbative calculations. Tabic 1 summa- 
rizes the fractional error for various physical quantities produced by 
truncating the e-expansion at leading, or next-to-leading, order before 
evaluating at e = 1, in the light Higgs limit. Although the lowest- 
order results often error by a factor of two or more, all but one of the 
next-to-leading order results are correct to better than 10%. (The free 
energy difference at the limit of metastability, AF(Tq), has the most 
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poorly behaved e-expansion. However, if one instead computes the log- 
arithm of this quantity, then the next-to-leading order result is correct 
to within 17%. The baryon violation rate is not shown because, due 
to the way its e-expansion was constructed, the result is trivially the 
same as the three-dimensional answer when Ai «C ^i- See ref. [6] for 
details.) 



observable ratio 


LO NLO 


asymmetric correlation length ^asym 
symmetric correlation length ^sym 
latent heat AQ 
surface tension a 
free energy difference AF{To) 


0.14 -0.06 
0.62 -0.08 
-0.23 0.04 
-0.40 -0.02 
-0.76 -0.44 



TABLE 1. The fractional error in the e-cxpansion results, 
when computing prefactors through leading order (LO) and 
next- to-leading order (NLO) in e, when Ai <C ffi. 



^ ^ 9^- When X/g'^ is 0(1), the three-dimensional loop expansion 
is no longer trustworthy. However, one may still test the stability 
of e-expansion predictions by comparing 0{e"') and 0(e"+^) predic- 
tions — provided, of course, one can evaluate at least two non-trivial 
orders in the e-expansion. For most physical quantities this is not 
(yet) possible; determining the lowest-order behavior of the prefactor 
in expansion (12) requires a one-loop calculation in the final effective 
theory together with a two-loop evaluation of the (solution to the) 
renormalization group equations. A consistent next-to-leading order 
calculation requires a two-loop calculation in the final theory together 
with three-loop renormalization group evolution. Althouth two-loop 
results for the effective potential and beta functions arc known, three 
loop renormalization group coefficients in the scalar sector are not 
currently available. Nevertheless, by taking suitable combinations of 
physical cpiantitics one can cancel the leading dependence on the scale 
ratio fi/T and thereby eliminate the dependence (at next-to-leading 
order) on the three loop beta functions. For example, the latent heat 
depends on the scale as AQ ~ (/Lt/r)^+^ while the scalar correlation 
length ^ ~ (/i/T)~^. Therefore, the combination ^^AQ cancels the 
leading fj,/T ~ e*^^^/^^ dependence and thus requires only two-loop in- 
formation for its next-to-leading order evaluation. The result of the 
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(rather tedious) calculation may be put in the form 



e^y^AQ = T'-'f{fl\,) l + 6 + 0{e 



(13) 



where S, the relative size of the next-to-leading order correction, is 
plotted in figure 4. The correction varies between roughly ±30% for 
(zero temperature) Higgs masses up to 150 GeV. This suggests that 
the e expansion is tolerably well behaved for these masses. For larger 
masses the correction does not grow indefinitely, but is bounded by 
80%, suggesting that the e expansion may remain qualitatively useful 
even when it does not work as well quantitatively. 




400 600 
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1000 



Figure 4- The relative size of the next-to-leading order correction to Casym^Q in the 
e-expansion. The values are given as a function of the (tree-level) zero-temperature Higgs 
mass in minimal SU(2) theory (iV = 2) with g = 0.63. 



4. Conclusions 

The most useful technique for studying the electroweak phase transition 
clearly depends on the range of value of the Higgs mass. For sufficiently 
light Higgs (Mh ^ 40 GeV) perturbation theory is adequate. At interme- 
diate masses (40 GeV < Mh < 80 GeV) numerical simulations have been 
quite effective. For heavier Higgs masses, the e-expansion appears promis- 
ing. 

Although considerable progress has been made toward a quantitative 
understanding of the electroweak phase transition, much remains to be 
done. With continuing efforts, numerical results will undoubtedly improve, 
particularly for Higgs masses above Mw. Calculations of additional physi- 
cal quantities at next-to-leading order in the e-expansion should definitely 
be performed to further confirm the reliability of the method. 

Nevertheless, it should be noted that quantitative understanding of the 
phase transition is already reasonably good in the range of Higgs masses 
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for which the transition is sufficiently discontinuous to be compatible with 
baryogenesis. Hence, the major source of uncertainty about the viability of 
electroweak baryogenesis appears to be the lack of knowledge about which 
theory is correct (extra Higgs, super symmetry, etc.) not the limitations in 
our ability to compute the thermodynamics of these theories. 
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